Airborne Electromagnetic Survey Research Articles

Island hydrogeology in the tropics: Constraining a 3D variable-density groundwater flow and solute transport model with geophysics

Saltwater intrusion is the greatest risk to coastal community water supplies where they are dependent on fresh groundwater as the main source of supply. For small, fractured bedrock island aquifers, the fresh groundwater lens dynamics and transition zone geometry are complex. This study investigated the impacts of projected increases in groundwater pumping on a fresh groundwater lens to evaluate changes to the lens geometry and localized up-coning from the deeper, more saline aquifers beneath a small bedrock island in the tropics. A combination of traditional hydrogeological datasets and an airborne electromagnetic survey were used to develop a three-dimensional density-dependent groundwater flow and solute transport model using the SEAWAT code. This investigation represents one of very few studies that have taken such an approach. The model was calibrated using observed groundwater hydraulic heads and chloride concentration data, and calculated chloride values based on bulk conductivity measurements determined from inverted geophysical data. A staged calibration approach was adopted, firstly assessing the time-average lens extent and geometry, and secondly considering the seasonal groundwater level response. In the calibration, the geophysical data helped constrain the lens geometry in the absence of hydraulic head and chloride data. The calibrated model was used to test scenarios where groundwater pumping rates were increased above the current demand of 452 m3/d, showing that the lens is likely to be stable, i.e., its available storage is not expected to contract excessively, for extraction rates of up to 3,000 m3/d. The combined use of geophysical data and a numerical modeling approach was advantageous in investigating the lens characteristics. It also demonstrated how these techniques can be used together to evaluate coastal water resources and to manage water supply risks for coastal communities. The study demonstrated that the freshwater lens can likely support the freshwater demands of the remote community, and is a preferred option compared with high-cost and more complicated options such as seawater desalination and managed aquifer recharge.

Three-Dimensional Separate and Joint Inversions of Multi-Component Frequency-Domain Airborne Electromagnetic Data: Synthetic Model Studies.

Airborne electromagnetic (AEM) surveys using airborne mobile platforms enable rapid and efficient exploration of areas where groundwork is difficult. They have been widely used in fields such as shallow resource exploration and environmental engineering. Three-dimensional AEM inversion is the main technique used in fine structural interpretation. However, most current methods focus on separate component data inversions, which limit the kinds of structures that can be recovered in the inversion results. To address this issue, a method for the robust 3D joint inversion of multicomponent frequency-domain AEM data was developed in this study. First, a finite element method based on unstructured tetrahedral grids was used to solve the forward problem of frequency-domain AEM data for both isotropic and anisotropic media. During inversion, a limited-memory quasi-Newton (L-BFGS) method was used to reduce the memory requirements and enable the joint inversion of large-scale multicomponent AEM data. The effectiveness of our algorithm was demonstrated using synthetic models for both isotropic and anisotropic cases, with 5% Gaussian noise added to the modeling data to simulate the measured data for separate and joint inversions. The results of the synthetic models show that joint inversion has advantages over separate inversion in that it enables the recovery of finer underground structures and provides a novel approach for the fine interpretation of frequency-domain AEM data.

Three-Dimensional Inversion of Induced Polarization Effects in Airborne Time Domain Electromagnetic Data Using the GEMTIP Model

This paper discusses the physical and mathematical principles of the airborne induced polarization (IP) method. The possibility of extracting information about the IP properties of rocks from airborne survey data has become a subject of active research recently. We introduce a method for the joint inversion of the airborne EM data into the electrical conductivity and IP parameters based on the generalized effective-medium theory of induced polarization (GEMTIP). We also present the results of the inversion of the airborne EM data collected over the Echum Project Area, in Northwestern Ontario, Canada, into 3D conductivity and chargeability models. Obtaining IP physical property models from an airborne geophysical survey may result in a paradigm change in mineral exploration by pulling more information and value from airborne EM surveys.

Fast Ice Thickness Distribution in the Western Ross Sea in Late Spring

AbstractWe present a 700 km airborne electromagnetic survey of late‐spring fast ice and sub‐ice platelet layer (SIPL) thickness distributions from McMurdo Sound to Cape Adare, providing a first‐time inventory of fast ice thickness close to its annual maximum. The overall mode of the consolidated ice (including snow) thickness was 1.9 m, less than its mean of 2.6 ± 1.0 m. Our survey was partitioned into level and rough ice, and SIPL thickness was estimated under level ice. Although level ice, with a mode of 2.0 m and mean of 2.0 ± 0.6 m, was prevalent, rough ice occupied 41% of the transect by length, 50% by volume, and had a mode of 3.3 m and mean of 3.2 ± 1.2 m. The thickest 10% of rough ice was almost 6 m on average, inclusive of a 2 km segment thicker than 8 m in Moubray Bay. The thickest ice occurred predominantly along the northwestern Ross Sea, due to compaction against the coast. The adjacent pack ice was thinner (by ∼1 m) than the first‐year fast ice. In Silverfish Bay, offshore Hells Gate Ice Shelf, New Harbor, and Granite Harbor, the SIPL transect volume was a significant fraction (0.30) of the consolidated ice volume. The thickest 10% of SIPLs averaged nearly 3 m thick, and near Hells Gate Ice Shelf the SIPL was almost 10 m thick, implying vigorous heat loss to the ocean (∼90 W m−2). We conclude that polynya‐induced ice deformation and interaction with continental ice influence fast ice thickness in the western Ross Sea.

A model of transmissivity and hydraulic conductivity from electrical resistivity distribution derived from airborne electromagnetic surveys of the Mississippi River Valley Alluvial Aquifer, Midwest USA

Groundwater-flow models require the spatial distribution of the hydraulic conductivity parameter. One approach to defining this spatial distribution in groundwater-flow model grids is to map the electrical resistivity distribution by airborne electromagnetic (AEM) survey and establish a petrophysical relation between mean resistivity calculated as a nonlinear function of the resistivity layering and thicknesses of the layers and aquifer transmissivity compiled from historical aquifer tests completed within the AEM survey area. The petrophysical relation is used to transform AEM resistivity to transmissivity and to hydraulic conductivity over areas where the saturated thickness of the aquifer is known. The US Geological Survey applied this approach to a gain better understanding of the aquifer properties of the Mississippi River Valley alluvial aquifer. Alluvial-aquifer transmissivity data, compiled from 160 historical aquifer tests in the Mississippi Alluvial Plain (MAP), were correlated to mean resistivity calculated from 16,816 line-kilometers (km) of inverted resistivity soundings produced from a frequency-domain AEM survey of 95,000 km2 of the MAP. Correlated data were used to define petrophysical relations between transmissivity and mean resistivity by omitting from the correlations the aquifer-test and AEM sounding data that were separated by distances greater than 1 km and manually calibrating the relation coefficients to slug-test data. The petrophysical relation yielding the minimum residual error between simulated and slug-test data was applied to 2,364 line-km of AEM soundings in the 1,000-km2 Shellmound (Mississippi) study area to calculate hydraulic property distributions of the alluvial aquifer for use in future groundwater-flow models.

Airborne Electromagnetic Survey over the Touro Copper VMS World Class Deposit (NW Spain): Geological and Geophysical Correlation

Electromagnetic (EM) methods belong to the main geophysical techniques used in the mineral exploration of massive sulphides. For selecting EM anomalies as possible massive sulphide targets, it is important to combine the geophysical results with other geological and/or geochemical techniques. In 2015, Atalaya Mining started a new mineral exploration project in the Touro Cu deposit, combining geological, geochemical (ore over 0.2% Cu), and geophysical techniques. The geophysical survey consisted of helicopter-borne EM using the versatile time-domain electromagnetic (VTEM™) max system operated by Geotech Ltd. with full-waveform processing. In total, 509 line-km of geophysical data were acquired during the survey that was completed in 2018. The results showed the massive sulphide Touro ore to be typically of the order of 0.25 ohm·m (4S/m conductivity) and host rock in the range of 1000–30,000 ohm·m, measured directly on the drill core. This modelling agreed well with the sub-horizontal dips observed for the known Touro ore bodies. The conductance modelled by the plate estimation of the VTEM data were also in good agreement with those provided by Geotech Ltd. and the resistivity/conductivity measurements we made on the massive sulphide samples from several Touro ore bodies. The combination of flat dips, good conductance, shallow depth, and, lastly, lack of conductive overburden or noneconomic conductive stratigraphy, i.e., graphitic shales and sulphide iron formation make the Touro project an ideal target for airborne electromagnetic prospecting. This paper presents the excellent correlation observed between the EM airborne anomalies and the massive sulphide blocks of the Touro copper deposit. Favourable factors contributing to the success of the survey were the high contrast in resistivity/conductivity between the massive sulphide Touro ore and the amphibolite host rock and minimal interference from “nuisance” conductors, such as graphitic shales.

Brief communication: The hidden labyrinth: deep groundwater in Wright Valley, Antarctica

Abstract. Since the 1960s, a deep groundwater system in Wright Valley, Antarctica, has been the hypothesized source of brines to hypersaline Don Juan Pond and Lake Vanda, both of which are rich in calcium and chloride. Modeling studies do not support other possible mechanisms, such as evaporative processes, that could have led to the current suite of ions present in both waterbodies. In 2011 and 2018, an airborne electromagnetic survey was flown over Wright Valley to map subsurface resistivity (down to 600 m) in exploration of liquid water. The surveys revealed widespread unfrozen brine in the subsurface near Lake Vanda, Don Juan Pond, and the North Fork of Wright Valley. While our geophysical survey can neither confirm nor deny deep groundwater connectivity between Lake Vanda and Don Juan Pond, it does point to the potential for deep valley-wide brine, likely within the Ferrar Dolerite formation.

Novel Airborne EM Image Appraisal Tool for Imperfect Forward Modeling

Full 3D inversion of time-domain Airborne ElectroMagnetic (AEM) data requires specialists’ expertise and a tremendous amount of computational resources, not readily available to everyone. Consequently, quasi-2D/3D inversion methods are prevailing, using a much faster but approximate (1D) forward model. We propose an appraisal tool that indicates zones in the inversion model that are not in agreement with the multidimensional data and therefore, should not be interpreted quantitatively. The image appraisal relies on multidimensional forward modeling to compute a so-called normalized gradient. Large values in that gradient indicate model parameters that do not fit the true multidimensionality of the observed data well and should not be interpreted quantitatively. An alternative approach is proposed to account for imperfect forward modeling, such that the appraisal tool is computationally inexpensive. The method is demonstrated on an AEM survey in a salinization context, revealing possible problematic zones in the estimated fresh–saltwater interface.

Motion noise remove for secondary-field semi-airborne transient electromagnetic data

Semi-airborne Transient Electromagnetics (SATEM) is a low-cost and productive method for quick geophysical surveys. The system includes a grounded line source on the ground and a receiver coil hanging in the air by a drone. Motion noise usually has the biggest impact on the acquired data because the receiver coil, which is suspended several meters below the drone by a soft rope, would maintain continuing random motion due to the movement of the drone, wind, etc. It also is a common problem in airborne electromagnetic (AEM) surveys, which typically use fixed-wing aircraft or helicopters. However, due to the secondary-field acquisition of SATEM, there are data deviations between the start and end of each motion period caused by continuous coil motion and discontinuous data acquisition. And this also brings difficulty in baseline elimination by using most de-noising methods, which require time or signal continuity on the baseline and are widely applied in full-time acquisition data (full-wave data). We proposed a simple but useful single-period polynomial-fitting method to remove motion noise from secondary-field data. And to prevent the obtained polynomial from being distorted from baseline by valid signals, a haar-based wavelet was employed to automatically recognize the start-point of pulses in the original signal so that the pulse can be eliminated before the polynomial-fitting process. Through a field SATEM survey data test, we found that the polynomial-fitting with 5 parameters (4th order) is good enough to describe the motion characteristics. We compared the results with those of the most commonly used wavelet methods at different bases and levels. The comparisons of de-noising results and apparent resistivity profiles of the original signal, data after wavelet de-nosing, and data after the polynomial-fitting de-noising show that this method is superior and much more reliable than wavelet de-noising of secondary-field SATEM data.

Improving spatial characterization of buried bedrock valleys through airborne frequency-domain electromagnetic, residual magnetic, and surface resistivity measurements

The development of a three-dimensional (3D) hydrogeologic site model requires detailed knowledge of an area's hydrostratigraphy, information typically obtained from core logs, historical well records, outcrop mapping, and/or surface geophysics. Due to the high cost of drilling, a 3D model may rely on various statistical techniques to interpolate lithological boundaries between points located tens of kilometers or more apart. This situation can result in poorly constrained or simplified lithostratigraphic models, with under-sampled features such as karst and buried bedrock valleys. The non-invasive nature of airborne geophysical surveys combined with its dense sampling capability can fill critical data gaps. This study evaluates the capacity of a helicopter borne frequency-domain electromagnetic (FDEM) and residual magnetic survey using the Resolve™ system over a buried bedrock valley incised in a multi-layered dolostone aquifer within a complex glaciated environment in an urban/agricultural region of southern Ontario, Canada. The resulting airborne data was evaluated against two high-resolution ground-based electrical resistivity transects oriented orthogonal to the valley, alongside a published Quaternary geological model. A comparison between our local geophysically-based model and the regional geological model shows that the airborne electromagnetic surveys can markedly improve characterization of bedrock valley morphology and Quaternary sedimentary unit architecture that is comparable to that obtained using high-resolution surface-based electrical resistivity measurements with 10 m electrode spacing. The scale of heterogeneity that can be resolved by the airborne survey reveals the level of petrophysical variability that can be seen within the valley infill. Our study demonstrates the potential of airborne FDEM surveys in semi-urban agricultural environments where detailed characterization of Quaternary sediment architecture and bedrock topography is needed to inform 3D groundwater flow system variability.

Joint estimation of groundwater salinity and hydrogeological parameters using variable-density groundwater flow, salt transport modelling and airborne electromagnetic surveys

Freshwater aquifers in low elevation coastal zones are known to be threatened by saltwater intrusion (SWI). As these areas host a significant share of the world's population, an excellent understanding of this phenomenon is required to effectively manage the availability of freshwater. SWI is a dynamic process, therefore saline groundwater distributions can change quickly over time – particularly in stressed areas with anthropogenic drivers. To model these changes, regional 3D variable-density groundwater (3D-VDG) flow and coupled salt transport models are often used to estimate the current (and future distributions) of saline groundwater. Unfortunately, parameterising 3D-VDG models is a challenging task with many uncertainties. Generally, uncertainty is reduced through the addition of observational data – such as Airborne Electromagnetic (AEM) surveys or ground-based information – that offer information about parameters such as salinity and hydraulic head. Recent research has shown the ability of AEM surveys to provide accurate 3D groundwater salinity models across regional scales, as well as highlighting the potential for good survey repeatability. To this end we investigated the novel approach of using repeat AEM surveys (flown over the same area at different points in time) and 3D-VDG models to jointly improve the parameterisation of 3D-VDG models - while simultaneously providing a detailed 3D map of groundwater salinity distributions. Using detailed 3D synthetic models, the results of this study quantitatively highlight the usefulness of this approach, while offering practical information on implementation and further research.

AEM in Norway: A Review of the Coverage, Applications and the State of Technology

From the first use of airborne electromagnetic (AEM) systems for remote sensing in the 1950s, AEM data acquisition, processing and inversion technology have rapidly developed. Once used extensively for mineral exploration in its early days, the technology is increasingly being applied in other industries alongside ground-based investigation techniques. This paper reviews the application of onshore AEM in Norway over the past decades. Norway’s rugged terrain and complex post-glacial sedimentary geology have contributed to the later adoption of AEM for widespread mapping compared to neighbouring Nordic countries. We illustrate AEM’s utility by using two detailed case studies, including time-domain and frequency domain AEM. In both cases, we combine AEM with other geophysical, geological and geotechnical drillings to enhance interpretation, including machine learning methods. The end results included bedrock surfaces predicted with an accuracy of 25% of depth, identification of hazardous quick clay deposits, and sedimentary basin mapping. These case studies illustrate that although today’s AEM systems do not have the resolution required for late-phase, detailed engineering design, AEM is a valuable tool for early-phase site investigations. Intrusive, ground-based methods are slower and more expensive, but when they are used to complement the weaknesses of AEM data, site investigations can become more efficient. With new developments of drone-borne (UAV) systems and increasing investment in AEM surveys, we see the potential for continued global adoption of this technology.

Current and future state of groundwater salinization of the northern Elbe-Weser region

Salinization of the upper aquifer of the northern Elbe-Weser region almost extends to the surface. Chloride content exceeds 250 mg/l and the groundwater is therefore, according to the German Drinking Water Ordinance, not suitable as drinking water. The chloride content in the aquifer originates from early flooding with seawater which occurred during the Holocene sea level rise. Depth and extent of the salinization were mapped by airborne electromagnetic surveys and validated by groundwater analyses. In the transition zone between the marshlands and geest areas, the fresh-saline groundwater interface falls to a depth of > −175 m NHN. Due to the extensive drainage of the marshlands, seepage of fresh groundwater is impeded. Instead, an upconing of the fresh-saline groundwater interface appears due to an upwardly directed hydraulic gradient. Due to climate change, chloride concentrations will increase along the coastlines. Further inland, a decrease of chloride content in near-surface groundwater will occur.

Nested anisotropic geostatistical gridding of airborne geophysical data

Airborne geophysical surveys routinely collect data along traverse lines at sample spacing distances that are two or more orders of magnitude less than between-line separations. Grids and maps interpolated from such surveys can produce aliasing; features that cross flight lines can exhibit boudinage or string-of-beads artifacts. Boudinage effects can be addressed by novel gridding methods. Following developments in geostatistics, a nonstationary nested anisotropic gridding scheme is proposed that accommodates sampling and data anisotropy in geophysical surveys. Computation is reduced by including anchor points throughout the interpolation region that contain localized anisotropy information that is propagated throughout the survey area with a smoothing kernel. Additional anisotropy can be required at certain locations in the region to be gridded. A model selection scheme is proposed that uses Laplace approximations for determining whether increased model complexity is supported by the surrounding data. The efficacy of the method is shown using a synthetic data set obtained from satellite imagery. A pseudogeophysical survey is created from the image and reconstructed with this method. Two case histories are selected for further elucidation from airborne geophysical surveys conducted in Western Australia (WA). The first example illustrates improvement in gridding the depth of paleochannels interpreted from along-line conductivity-depth models of a regional airborne electromagnetic (AEM) survey in the Mid West WA region. The second example indicates how improvements can be made in producing grids of aeromagnetic data and inverted electrical conductivity from an AEM survey conducted in the Pilbara region. In both case histories, nested anisotropic kriging reduces the expression of boudinage patterns and sharpens crossline features in the final gridded products permitting increased confidence in interpretations based on such products.

1D Stochastic Inversion of Airborne Time-Domain Electromagnetic Data with Realistic Prior and Accounting for the Forward Modeling Error

Airborne electromagnetic surveys may consist of hundreds of thousands of soundings. In most cases, this makes 3D inversions unfeasible even when the subsurface is characterized by a high level of heterogeneity. Instead, approaches based on 1D forwards are routinely used because of their computational efficiency. However, it is relatively easy to fit 3D responses with 1D forward modelling and retrieve apparently well-resolved conductivity models. However, those detailed features may simply be caused by fitting the modelling error connected to the approximate forward. In addition, it is, in practice, difficult to identify this kind of artifacts as the modeling error is correlated. The present study demonstrates how to assess the modelling error introduced by the 1D approximation and how to include this additional piece of information into a probabilistic inversion. Not surprisingly, it turns out that this simple modification provides not only much better reconstructions of the targets but, maybe, more importantly, guarantees a correct estimation of the corresponding reliability.

Thermal legacy of a large paleolake in Taylor Valley, East Antarctica, as evidenced by an airborne electromagnetic survey

Abstract. Previous studies of the lakes of the McMurdo Dry Valleys have attempted to constrain lake level history, and results suggest the lakes have undergone hundreds of meters of lake level change within the last 20 000 years. Past studies have utilized the interpretation of geologic deposits, lake chemistry, and ice sheet history to deduce lake level history; however a substantial amount of disagreement remains between the findings, indicating a need for further investigation using new techniques. This study utilizes a regional airborne resistivity survey to provide novel insight into the paleohydrology of the region. Mean resistivity maps revealed an extensive brine beneath the Lake Fryxell basin, which is interpreted as a legacy groundwater signal from higher lake levels in the past. Resistivity data suggest that active permafrost formation has been ongoing since the onset of lake drainage and that as recently as 1500–4000 years BP, lake levels were over 60 m higher than present. This coincides with a warmer-than-modern paleoclimate throughout the Holocene inferred by the nearby Taylor Dome ice core record. Our results indicate Mid to Late Holocene lake level high stands, which runs counter to previous research finding a colder and drier era with little hydrologic activity throughout the last 5000 years.

Stratigraphic identification with airborne electromagnetic methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell's equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

Induced polarization effects in airborne transient electromagnetic data collected in the McMurdo Dry Valleys, Antarctica

SUMMARYAirborne electromagnetics (EM) is a geophysical tool well suited to mapping glacial and hydrogeological structures in polar environments. This non-invasive method offers significant spatial coverage without requiring access to the ground surface, enabling the mapping of geological units to hundreds of metres depth over highly varied terrain. This method shows great potential for large-scale surveys in polar environments, as common targets such as permafrost, ice and brine-rich groundwater systems in these settings can be easily differentiated because of their significant contrasts in electrical properties. This potential was highlighted in a 2011 airborne EM survey in the McMurdo Dry Valleys that mapped the existence of a large-scale regional groundwater system in Taylor Valley. A more comprehensive airborne EM survey was flown in November 2018 to broadly map potential groundwater systems throughout the region. Data collected in this survey displayed significant perturbations from a process called induced polarization (IP), an effect that can greatly limit or prevent traditional EM workflows from producing reliable geological interpretations. Here, we present several examples of observed IP signatures over a range of conditions and detail how workflows explicitly designed to handle IP effects can produce reliable geological interpretations and data fits in these situations. Future polar EM surveys can be expected to encounter strong IP effects given the likely presence of geological materials (e.g. ice and permafrost) that can accentuate the influence of IP.

Linking sea ice deformation to ice thickness redistribution using high-resolution satellite and airborne observations

Abstract. An unusual, large, latent-heat polynya opened and then closed by freezing and convergence north of Greenland's coast in late winter 2018. The closing presented a natural but well-constrained full-scale ice deformation experiment. We observed the closing of and deformation within the polynya with satellite synthetic-aperture radar (SAR) imagery and measured the accumulated effects of dynamic and thermodynamic ice growth with an airborne electromagnetic (AEM) ice thickness survey 1 month after the closing began. During that time, strong ice convergence decreased the area of the refrozen polynya by a factor of 2.5. The AEM survey showed mean and modal thicknesses of the 1-month-old ice of 1.96 ± 1.5 m and 1.1 m, respectively. We show that this is in close agreement with modeled thermodynamic growth and with the dynamic thickening expected from the polynya area decrease during that time. We found significant differences in the shapes of ice thickness distributions (ITDs) in different regions of the refrozen polynya. These closely corresponded to different deformation histories of the surveyed ice that we reconstructed from Lagrangian ice drift trajectories in reverse chronological order. We constructed the ice drift trajectories from regularly gridded, high-resolution drift fields calculated from SAR imagery and extracted deformation derived from the drift fields along the trajectories. Results show a linear proportionality between convergence and thickness change that agrees well with the ice thickness redistribution theory. We found a proportionality between the e folding of the ITDs' tails and the total deformation experienced by the ice. Lastly, we developed a simple, volume-conserving model to derive dynamic ice thickness change from the combination of Lagrangian trajectories and high-resolution SAR drift and deformation fields. The model has a spatial resolution of 1.4 km and reconstructs thickness profiles in reasonable agreement with the AEM observations. The modeled ITD resembles the main characteristics of the observed ITD, including mode, e folding, and full width at half maximum. Thus, we demonstrate that high-resolution SAR deformation observations are capable of producing realistic ice thickness distributions.

Suppression of motion noise based on a linear‐homomorphic filtering algorithm in airborne electromagnetic survey

ABSTRACTThe movement of the receiver coil in a magnetic field produces motion noise in the airborne electromagnetic field. With an amplitude that is several orders of magnitude higher than that of other noise sources, this motion noise is one of the main low‐frequency (<1 kHz) noise sources. Especially in the late decay period, the residual motion noise produces false anomalous information, which is a key problem of low‐base‐frequency airborne electromagnetic surveys. In this paper, a method for calculating the motion noise is proposed. This method takes into account both the earth's magnetic field and the active secondary magnetic field and introduces additive and multiplicative motion noise. Then, a linear‐homomorphic filtering algorithm is proposed to handle the additive and multiplicative motion noises. In the linear filtering part of our algorithm, the minimum noise fraction transform is applied to suppress additive motion noise. It provides sets of components ordered by the signal‐to‐noise ratio by deriving the eigenvalues and eigenvectors for the signal and uncorrelated noise covariance matrices. In the homomorphic filtering part, the multiplicative motion noise is removed. To do this, the adaptive‐width smoothing algorithm is applied to the logarithm of the filtered with the linear filtering process profile data. The results of synthetic and field data demonstrate that the linear filtering process based on the minimum noise fraction transform is effective for very low signal‐to‐noise ratios of the late‐time data, and the homomorphic filtering process can help distinguish small anomalies of these channels.